Interplay between writhe and knotting for swollen and compact polymers 
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The role of the topology and its relation with the geometry of biopolymers under different phys- 
ical conditions is a nontrivial and interesting problem. Aiming at understanding this issue for a 
related simpler system, we use Monte Carlo methods to investigate the interplay between writhe 
and knotting of ring polymers in good and poor solvents. The model that we consider is interacting 
self-avoiding polygons on the simple cubic lattice. For polygons with fixed knot type we find a 
writhe distribution whose average depends on the knot type but is insensitive to the length iV of the 
polygon and to solvent conditions. This "topological contribution" to the writhe distribution has 
a value that is consistent with that of ideal knots. The standard deviation of the writhe increases 
approximately as in both regimes and this constitutes a geometrical contribution to the writhe. 
If the sum over all knot types is considered, the scaling of the standard deviation changes, for com- 
pact polygons, to ~ N ' 6 . We argue that this difference between the two regimes can be ascribed 
to the topological contribution to the writhe that, for compact chains, overwhelms the geometrical 
one thanks to the presence of a large population of complex knots at relatively small values of N. 
For polygons with fixed writhe we find that the knot distribution depends on the chosen writhe, 
with the occurrence of achiral knots being considerably suppressed for large writhe. In general, the 
occurrence of a given knot thus depends on a nontrivial interplay between writhe, chain length, and 
solvent conditions. 



PACS numbers: 02.10.Kn, 36.20.Ey, 87.15.A-, 36.20.- 



I. INTRODUCTION 

Important biopolymers such as duplex DNA exist as 
double-stranded macromolecules, where the two strands 
of complementary nucleotides are wound around each 
other in a right-handed fashion and around a common 
axis In addition, the double helix can wind in space 
to form a new super-helix, in which case the polymer is 
said to be supercoiled. A geometric quantity that has 
proved to be useful in describing the degree of supercoil- 
ing in double stranded DNA (dsDNA) is the writhe of 
a curve, and there is an important conservation theorem 
0, [|[ relating the writhe of the central axis of circular ds- 
DNA to the double helical twist and the linking number 
of the two strands. Supercoiling in DNA can result from 
its binding to proteins (histones) in chromatins or from a 
linking deficit between the two strands, when the macro- 
molecule becomes a ring, for instance by cyclization. For 
example, circular DNA extracted from cells has negative 
supercoiling [ij. On the other hand random cyclization 
of linear dsDNA with cohesive ends can trap DNA topoi- 
somers as knots [E @. Both supercoiling and knotting 
in circular DNA are critical to the functioning of the cell 
and, for this reason, there exist enzymes to control 
the writhe and the knotting of DNA, especially during 
replication, transcription and recombination. This sug- 
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gests that supercoiling and knot formation are relevant 
indicators of the spatial arrangement of circular DNA 
and that a detailed study of their reciprocal influence is 
crucial in understanding the conformational properties of 
this macromolecule at equilibrium [E [|| . 

A relationship between writhe and knot type is evident 
also in the knot theory of random closed curves where it 
is known for example that chiral knots have a non neg- 
ligible writhe whose sign depends on the chirality of the 
knot, and that there are families of knots (e.g. torus and 
twist knots) that can be clearly distinguished in terms 
of their average writhe EE EH ■ For physical knots this 
relationship is more stringent and allows a classification 
of knots in terms of the geometrical (writhe) properties 
of their ideal representations EE EE EH . 

The above mentioned examples refer in general to sit- 
uations in which the polymer is in a swollen phase due 
to good solvent or unconstrained conditions. There are 
however several important cases in which macromolecules 
such as DNA and proteins are in highly condensed phases 
due to bad solvent conditions or to strong confinement. 
Highly compact configurations have geometrical proper- 
ties that clearly differ from their swollen counterparts 
and this may change the relationship between writhe and 
knotting quite dramatically. This is the case, for exam- 
ple, for highly condensed DNA extracted from P4 phages 
El EE EE EE d ■ The genome of the P4 phage is a lin- 
ear dsDNA, about 11 Kbp l ong and having two 19bp 
single stranded cohesive ends [19| . In P4 tailless mutants 
the two ends can move freely within the capsid and at 
the end of the packing process they can cohere to form 
rings. Since the genome is very long compared to the 
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linear size of the capsid (50nm) the DNA is very con- 
densed and a high probability of knotted molecules is 
expected. The experimental data show that over 97% 
of the molecules in tailless mutants are indeed knotted. 
This is dramatically higher than the 3% value observed 
when P4 DNA molecules undergo cyclization in free so- 
lution [f| [2(J. Moreover, knots in P4 phages are very 
complex, being characterized by a crossing number of- 
ten much greater than 10 [13] • A more detailed analysis 
of the extracted knot spectrum based on 2D gel elec- 
trophoresis reveals a predominance of chiral knots over 
achiral ones and a prevalence of torus knots over twist 
knots [H, [2l[ . These findings suggest a spatial organiza- 
tion of the viral genome characterized by a large amount 
of writhe that induces a bias in the observed knot spec- 
trum towards chiral and, especially, torus knots fl8l. |2l|. 
This may suggest an important contribution of the writhe 
to the condensation of DNA [22j • 

All the examples outlined above suggest the impor- 
tance of a theoretical study of the interplay between 
writhe and knot distribution for polymer rings under dif- 
ferent equilibrium conditions. Here we explore the re- 
lation between writhe and knotting for the two extreme 
cases of polymers in the swollen and compact phases. 

We perform Monte Carlo simulations on a simple dis- 
crete model of a self-attracting ring polymer, namely 
a polygon in the cubic lattice in which the quality of 
the solvent is mimicked by assigning an energy gain to 
any pair of nearest neighbor non-consecutive vertices of 
the polygon. By tuning properly the energy gain this 
model is known to describe a collapse (8) transition 
from a swollen phase, characterized by extended con- 
figurations, to one in which the polygons are compact 
[23,[24, 25]. This allows us to compare, within the same 
model, the relationship between writhe and knotting in 
the two regimes and to see, for example, to what ex- 
tent the degree of condensation of the polymer affects 
the writhe distribution and its relation to the knot spec- 
trum. 

The plan of this paper is as follows. In the next section 
we describe the model, the Monte Carlo algorithms and 
the techniques adopted to compute the writhe and to 
identify the knot type. In Sec. IIIII we present the results 
first for the writhe distribution at fixed knot type and 
then for the knot spectrum as a function of the writhe. 
Section lTVl is devoted to a general discussion of the results 
presented and to some conclusions. 



II. MODEL AND SIMULATION METHODS 

A. Polymer model 

As a model for the large scale behavior of a long flexible 
circular polymer in a good solvent (swollen phase), we use 
TV-step self avoiding polygons (SAPs) on the cubic lattice, 
i.e. closed lattice walks whose steps can visit each vertex 
of the lattice at most once [23[ . To mimic the quality of 



the solvent we add to the model an effective attractive in- 
teraction potential which lowers the total energy by e = 1 
whenever two non-consecutive vertices of the SAP are 
one lattice distance apart. This attractive interaction is 
sufficient to induce a collapse transition at a critical tem- 
perature T c and, for T < T c , it is known to describe the 
compact phase [231 ]. In this work we have used T = 2.5, 
well below the value T c ~ 3.717 [H, HH]. A collapsed 
configuration with TV = 1000 steps is shown in Fig.QJa). 
Equilibrium configurations are sampled by Monte Carlo 
simulations based on two different algorithms, one for 
each equilibrium phase. In the swollen phase the algo- 
rithm uses the set of two-pivots moves i.e. non-local defor- 
mations that are known to sample efficiently unweighted 
polygons in Z 3 [26j ] . 

In the compact phase the two-pivots algorithm is 
known to be inefficient [25| and in order to have good 
statistics in this regime we use the pruned enriched 
Rosenbluth method (PERM) . This is a walk growing al- 
gorithm that has been shown to be quite effective in sam- 
pling compact self-avoiding walks [27| . Note that here the 
efficiency of the algorithm is mitigated by the fact that 
only closed chains, i.e. a small fraction of the whole set of 
linear chains, are accepted. Nevertheless this limitation 
is not very severe in the compact phase and we have been 
able to sample compact polygons with N ranging from 
200 (680000 configurations) up to 1600 (220000 configu- 
rations) . For this range of N the compact phase is char- 
acterized by a sufficiently high knotting probability and 
knot complexity [28|. This is not the case in the swollen 
phase where the knot population becomes non negligible 
only for large values of TV. In this case, in order to have 
good statistics of knots, we sample configurations with 
N ranging from 1000 up to 200000. 



B. Computation of the writhe and identification of 

knots 

A commonly used algorithm to compute the writhe of 
a curve goes as follows: First, one projects the curve 
onto an arbitrary plane. In general the projection will 
have crossings that most of the time will be transverse, 
so that, after having established an orientation of the 
curve, a sign +1 or —1 (determined by a right hand rule) 
can be assigned to each crossing. The sum of these signs 
gives the signed crossing number in this projection. The 
writhe of the curve is obtained by averaging these signed 
crossing numbers over all possible projections. From this 
definition it is clear that the main difficulty in comput- 
ing the writhe of a configuration would be the averaging 
procedure over all projections. Fortunately, for polygons 
in Z 3 this procedure is enormously simplified by a the- 
orem [29J which reduces the writhe computation to the 
average of linking numbers of the given curve with four 
selected push-offs of the curve itself. In our calculation 
of the writhe we made extensive use of this result. 

The determination of the knot type of a given config- 
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FIG. 1: (a) Example of a N = 1000 steps compact configuration and its projection on a plane, (b) The configuration of 
Fig. QJa) after the simplification procedure based on the BFACF algorithm. Note that in this case the simplification has been 
able to reach the minimal number of steps (N = 24) compatible with the knot type (3i) of the polygon (41J. 



uration, is in general, a difficult task both in the swollen 
phase where N is very large and in the compact phase 
where polygons are highly condensed. In both cases any 
projection of the chain onto a plane gives rise to a knot 
diagram with many crossings that makes the calculation 
of polynomial invariants quite prohibitive [see Fig.fija)]. 
To circumvent this difficulty, we simplify each sampled 
configuration before performing the projection. This is 
achieved by applying to the polygon a smoothing algo- 
rithm which progressively reduces the length of the chain 
while keeping its knot type unaltered (for a similar pro- 
cedure, see [308) ■ This procedure is based on the BFACF 
algorithm [3l|, [32|, an A*"- varying Monte Carlo method 
that is known to be ergodic within each knot type and 
in which the step fugacity, if kept low enough, induces a 
rapid reduction in the number of edges in the polygon. 
This simplification technique can reduce dramatically the 
number of crossings encountered in an arbitrary projec- 
tion, as shown in Fig. [IJb) for the case of Fig. [IJa) and 
more generally in Fig.[2]for compact and swollen polygons 
with A" = 1000. Note that the simplification procedure 
has a dramatic effect on the number of crossings. 

For each simplified configuration we perform 500 pro- 
jections and we choose the projection with the minimal 
number of crossings. The resulting knot diagram is en- 
coded in terms of the Dowker code [33j . A further sim- 
plification of the Dowker code based on Reidemeister- 
like moves is performed. Finally, a factorization of the 
Dowker code is attempted. This procedure, whenever 
successful, splits composite knots into their prime com- 
ponents. From each component of the original Dowker 
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FIG. 2: Histogram (with arbitrary normalization) of the av- 
eraged (over 500 projections) number of crossings before and 
after the simplification procedure. Panels (a) and (b) refer 
to polygons with N = 1000 respectively in the compact and 
swollen phase. 



code we extract, by using Knotfind [34[, the knot type 
of the original configuration. In this way we have been 
able to distinguish composite knots with up to 5 prime 
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FIG. 3: The dense histogram is the distribution of the writhe 
after the simplification based on the BFACF algorithm for 
the knot type 3i (N = 1600, T = 2.5). This is divided in two 
sectors Wr < and Wr > to emphasize the detection of the 
two chiralities. Distributions of the writhe corresponding to 
original configurations for each chirality are shown as empty 
histograms. Bins have size 1/4, i.e. the minimum resolution 
of Wr on the cubic lattice. 



components, and with each component having crossing 
number up to 11. 

To identify the chirality of a given knot type we can- 
not simply rely on the above procedure, since the Dowker 
code does not account for the handedness of the knot. 
Following [35| we used a heuristic approach based on the 
well known correlation between the chirality of a knot 
and the writhe of its minimal diagrammatic representa- 
tion [T3, H3 • This corresponds to computing the distri- 
bution of the writhe of the simplified configurations and 
looking at the shape of this distribution. If the distri- 
bution is sharply peaked around a well defined positive 
or negative value of the writhe one can assert with confi- 
dence that the knot is chiral with corresponding sign. In 
Fig. [3] we show the distribution of the writhe for the sim- 
plified configurations with a knot 3i, for N — 1600 and 
T = 2.5. The distribution is clearly bimodal, and one can 
reliably distinguish most of the knot chiralities. Distribu- 
tions of corresponding configurations before the BFACF 
procedure are also shown for comparison. Clearly the 
main uncertainties in this method of chiral detection are 
due to configurations whose writhe, after BFACF simpli- 
fication, has a value close to zero. For these, relatively 
few, configurations we have calculated the Jones polyno- 
mial, and hence their chirality, explicitly. 



III. NUMERICAL RESULTS 

A. The writhe distribution for fixed knot type 

We first focus on the writhe distribution for config- 
urations with fixed knot type r, Ppf(Wr\r). In Fig. [3] 
and Fig.[5]we show the distributions Pjv( Wr|r) for some 
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FIG. 4: Distributions P(Wr\r) in the compact phase for 
N = 200 (narrower), 400, 800, and 1600: (a) r = 3f , (b) 3J", 
(c) 4i, (d) 5+, and (e) 5^. 



knots respectively in the compact and swollen phases. In 
each panel different curves correspond to different values 
of N. We stress again that, in the good solvent regime, 
the knot probability is very low unless N is large. Conse- 
quently we have to consider polygons with N values that 
are up to two orders of magnitude larger than their com- 
pact counterparts. For the 3i knot both mirror images 
are shown. The plots suggest a very small (if any) de- 
pendence on N for the average writhe /i(r) whereas as N 
increases the distributions tend to broaden. This behav- 
ior, already observed in [HI HH for swollen polygons is 
here confirmed also for compact configurations. The dif- 
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FIG. 5: Distributions P( Wr\r) in the swollen phase for N = 
5000 (narrower), 10000, 30000, and 80000: (a) unknot (b) 
r = 3+ (c)3r/, (d)4i. 



ferent ranges in Wr between swollen and compact phases 
are due to the different values of N considered in the two 
situations. If the same value of N is considered (see Fig. [6] 
for N — 1000) the distribution P^(Wr\r) looks similar 
in the two phases, being slightly broader for compact 
polygons. 

In Fig. [7] we show the writhe distribution for achiral 
compact polygons for N — 1600: all the curves are sym- 
metric and centered around zero. It was shown rigorously 
in [Icl ] that achiral knots have mean writhe equal to zero 
for un- weighted polygons. That argument is a symmetry 
argument that can be extended mutatis mutandis to the 
case of compact polygons. 
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FIG. 6: Distributions P( Wr\r) for polygons with N = 1000, 
both in the swollen (thick dashed line) and in the compact 
(dashed line) phase . The thin solid lines correspond to one 
parameter Gaussian fits in which the value of the average is 
fixed to the ideal value /i(3*) ~ 3.41. 
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FIG. 7: Distributions P( Wr\r) for polygons with N = 1600 
in the compact phase, for achiral knots r = 4i and r = 63, 
and for a composite achiral knot r = 3^#3f. The solid line 
is a one parameter Gaussian fit with zero mean. 



In Fig. [8] we report the N dependence of the average 
writhe for different knot types, respectively in the com- 
pact and swollen phase. Within error bars (68% con- 
fidence intervals) it is clear that the average writhe is 
essentially independent of N and the values are in good 
agreement with those calculated previously for different 
models of polymers in good solvent [H, Q and for ideal 
knots [llilrj] (dashed lines in the plots). The agreement 
between the average writhe for a fixed knot type r and 
the value for the corresponding ideal knot is extended to 
composite knots as shown in Fig. [5] indicating that the 
additivity property of the writhe under knot composition 
found for ideal knots in [36[ and for unweighted polygons 
in [l(| applies also for the case of compact polygons. In 
fact Sumners [37| has recently proved a result about ad- 
ditivity of writhe under the connect sum operation for 
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FIG. 9: Mean writhe for knots 3i and all the possible composite knots 3i#3i in the compact (left panel) and swollen (right 
panel) phase. Dashed lines are the values of the sum of the means writhe of the corresponding ideal prime knots. 



unweighted polygons. 

The fact that the average writhe is independent of the 
quality of the solvent is a further indication that this 
number reflects topological properties of the rings no 
matter how badly they are embedded in space. One can 
understand this property by the following argument: for 
a given configuration and projection one may attempt to 
simplify the resulting knot diagram by performing Rei- 
demeister moves with the goal of removing most of the 
inessential crossings. This will give a diagram with the 
minimal number of crossings compatible with the projec- 
tion and with the given knot type r. Note that in the 
cubic lattice this procedure would correspond to a sim- 
plification of the original configuration for example by 
using the BFACF algorithm. This will bring the original 
configuration, either swollen or compact, close to the one 
with the minimal number of steps (ideal configuration 
on the cubic lattice). Let us now see how this procedure 
affects the average writhe h{t). Since the Reidemeis- 
ter II and Reidemeister III moves always involve pairs of 
crossings with opposite signs, the planar writhe is unaf- 



fected. This is not true for the first Reidemeister move 
that removes single crossings. For a particular conforma- 
tion the number of crossings with positive and negative 
signs that can be removed by Reidemeister I moves will 
not be equal. However, the procedure of averaging over 
all conformations should drastically reduce this average 
difference. The resulting average planar writhe is then 
related to the essential crossings for a given knot type 
and should coincide with the value for the ideal knot. 
Note that the above arguments suggest that the distri- 
bution P/v(WV|t) has Gaussian tails. In fact if the knot 
is achiral (as in Fig. [7]) or if chiral knots are split into mir- 
ror components (as for example the plus trefoil in Fig.O, 
the writhe distribution turns out to be well approximated 
by a Normal distribution Af (/x(t), <t(t)) with mean /Lt(r) 
and standard deviation a(r) (see solid lines in Fig. [6] and 
Fig.©. 

By looking at the width of the distributions in Fig- 
ures H] and [5] it is clear that the standard deviation (Jn(t) 
increases with N and this dependence is shown in Fig. [Tffl 
left panel, for compact prime knots and for swollen knots 
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FIG. 10: Log-log plot of the standard deviation <tjv(t) of the writhe distribution, as a function of N, for several knots (averaged 
over mirror images), in the compact (left panel) and swollen (right panel) phase. For each knot type the data have been fit 
by a power law of the form A(t)N Vt (dashed lines) with r\ T = 0.525 ± 0.020 in the compact phase and r\ T — 0.50 ± 0.01 in 
the swollen one. Note that for compact polygons, in order to avoid correction to scaling effects, the fit has been obtained by 
excluding the data N = 200, 400. 



in the right panel. A simple fit of the data in the form 
A(t)N Vt gives estimates of the exponents that, within er- 
ror bars, are independent on the knot type r and whose 
value is 0.525±0.020 for compact polygons and 0.50±0.01 
for swollen ones. Note that the two values are, within 
error bars, compatible suggesting that r\ T is also inde- 
pendent of the equilibrium phase, and possibly equal to 
1/2- 

Given that r\ r seems to be independent of r it is inter- 
esting to see if this is also the case for the amplitude A(t). 
In Fig.[TT]we plot <jn{t)/N 11t for different knot types, re- 
spectively in the compact (left panel) and swollen (right 
panel) phase. Within error bars the amplitude turns out 
to be independent of the knot type. 

It is interesting to compare the TV dependence of 
itjv(t) [i.e. the standard deviation of the distributions 
P/v(W / r|r)] with the one, cat, that measures the broad- 
ness of Pn{Wt) = J2 T ^n(t)Pn( Wt|t) where 7Tat(t) is 
the probability of occurrence of knot type r. Note that 
Pn ( Wr) is the writhe distributions of the whole set of 
polygons. In Fig. [T^] we plot ctjv as a function of N 
for swollen (empty circles) and compact (empty squares) 
polygons. Moreover, since often in the literature, the 
width of the writhe distribution Wr), is described in 
terms of the mean of the absolute value of the writhe 
(I VFr|)jv, in Fig. [12] we have also reported this quantity 
in terms of solid symbols. The plots are on a log-log scale 
and the linear behavior suggests a power law dependence 
AN V . Note that for each phase solid and empty sym- 
bols lie on parallel lines indicating that ctjv and (| Wr\) jv 
display the same power law behavior but different am- 
plitude. A simple linear fit of the data furnishes the 
estimate r\ = 0.5035 ± 0.0006 in the swollen phase and 
1] = 0.599 ± 0.002 in the compact one. The value found 
for the swollen regime agrees with previous estimates [38[ 
and is close to the rigorous lower bound 1/2 established 
for (| Wr\) 38]. In the compact phase, however, the expo- 



nent r\ is significantly higher than its swollen counterpart 
but mark edly smaller than the rigorous upper bound 4/3 
found in [39|. Its value is also smaller than 0.75, i.e. the 
one estimated in [3(3], for Gaussian chains confined in 
small spheres. The difference between these two last es- 
timates could be due to the presence, in our case, of 
the excluded volume interactions that increase the per- 
sistence length of the chain making the wrapping of the 
chain around itself more costly. However, the difference 
in the value of 77 could be due to a difference in the spa- 
tial organization of compact configurations obtained by 
self- attraction compared to the ones obtained by a severe 
geometrical confinement. 

An important byproduct of the analysis of ctjv and 
<tjv(t) is that r\ w r\ T in the swollen phase whereas 77 > r\ T 
in the compact one. This should be contrasted with the 
fact that 7] T is nearly the same ~ 1/2 in the two phases. 

A possible explanation of this striking behavior may be 
that the writhe of a given configuration can be decom- 
posed into a "geometrical" contribution coming from the 
spatial organization of the polygon and the topological 
one coming purely from the knot type of the configura- 
tion. For the swollen case it is known that the geometri- 
cal contribution to a increases at least as rapidly as yN 
which is what was found for the set of all polygons in 
the swollen phase. In this phase, the topological con- 
tribution apparently does not scale faster than y/N. In 
general however the situation could be different since a 
sum over all knot types must be taken into account and 
the topological contribution could become important if 
the probability of a knot type r to occur, i.e. ttn(t) is 
sufficiently large. This turns out to be the case for com- 
pact polygons where the knotting probability increases 
with N with an exponential factor that is a factor of 
1000 bigger than for its swollen counterparts. 

To make these arguments more precise let us suppose 
that P/v(WV|t) are well approximated by Gaussian dis- 
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FIG. 11: Standard deviation <t(t) divided by N Vt for different knot types in the compact (left panel) and swollen (right panel) 
phase. The value of rj T considered is 0.525 ± 0.020 for compact polygons and 0.50 ± 0.01 for the swollen ones. 



tributions. If this is the case (as suggested by our re- 
sults) the conditional distributions PN(Wr\r) are com- 
pletely defined by the average ^n(t) and the standard 
deviation cjn{t)- If knots are distinguished by chirality, 
clearly the variance of the distribution is £t/v(t) 2 . If we 
sum over the mirror images the mean writhe is zero and 
the overlap of two Gaussian distributions implies that 
(Wr 2 )N{T) = (Jn(t) 2 + hn{t) 2 - This term is for fixed 
knot type and, in order to compute the global second 
moment, the term must be weighted by the probability 
of occurrence of the knot 

(Wr 2 ) =5> jv (t)(WV 2 ) jv (t). 

T 

This gives 

(Wr 2 )^(a N (r) 2 ) + ( m (r) 2 ). 

The scaling of (Wr 2 ) depends on whether (a 2 ) or (/1 2 ) 
increases more rapidly as N increases. For the compact 
case it seems that (/1 2 ) increases more rapidly than (a 2 ) 
while this is not true in the swollen phase. We can es- 
timate (/iAr(r) 2 ) in the swollen phase by using the ideal 
Hn(t) whenever a knot r contributes to the statistics. 
Since the two members of a chiral pair occur with roughly 
equal probabilities and since knots are weakly localized in 
the swollen phase we can estimate the writhe of a compos- 
ite knot (for all chiral combinations) by assuming inde- 
pendence. The average ((^n ( t ) 2 )) 1 ^ 2 from this analysis 
turns out to scale ~ V~N, see Fig. [TJ1 and with an ampli- 
tude much smaller than cr/v, confirming that the scaling 
of ctjv is governed by the geometrical spread (o-jv(t) 2 ) of 
all knots around their ideal value. 

Hence, the important difference between the swollen 
phase and the compact one is that collapsed polymers 
display a rich spectrum of knots n at(t) for moderate 
lengths N 28], while in good solvent conditions knots 
start to appear with a non-negligible frequency only for 
N > 10 5 [4(3 . At that stage cr(r) > |//(r)| and hence 
(M / r 2 )Ar(r) ~ er(-r) 2 . On the contrary, in the compact 
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FIG. 12: Log-log plot of u(t) vs in the swollen phase (cir- 
cles) and in the compact phase (squares). The filled sym- 
bols refer to the mean of the absolute value of the writhe. 
Dashed lines are the result of a linear ht and have a slope of 
0.599 ± 0.002 and 0.5035 ± 0.0006 respectively for compact 
and swollen polygons. Stars refer to the average square of the 
mean value of writhe of knots in the swollen phase and the 
dotted line is a ht with slope 0.499 ± 0.002. 



phase knots start to appear with a non-negligible fre- 
quency already for A^'s where |/x(r)| > ^(t), in which 
case (Wr 2 )N(r) ~ m( t ) 2 is a good approximation. Be- 
cause of this effect, the global (Wr 2 ) has a scaling that 
does not trivially follow that of single knots. 



B. Knot spectrum at fixed writhe 

Podtelezhnikov et al [8| have investigated the prob- 
ability of finding a specific knot type, conditioned on 
the writhe, for a model of circular DNA in a good 
solvent. In this section we examine how the pres- 
ence of writhe might affect the knot spectrum in 
the swollen and compact phases. The idea is to 
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extract a subset of the sample of polygons whose 
writhe falls into a window of values [WV m ; n , Wr max ]. 
Within this subset one looks at the relative frequency 
of occurrence of a given knot. This corresponds 
to estimating the probability P(r|[W r r min , Wr max ]) = 

C:::: m w r y e t j^t;; Wr )- . In Fi s- we 

plot the knot spectrum for two different intervals of the 
writhe (middle and lower panels) for polygons in the com- 
pact phase. The top panel is the knot spectrum obtained 
by summing over all values of the writhe and is shown as 
a reference. Let us focus first on the N = 600 case (empty 
bars). Here it is clear that, as writhe increases, the knot 
spectrum changes dramatically with respect to the case 
of unconstrained writhe (top panel). In particular one 
sees that in configurations with relatively large value of 
writhe (bottom panel) the achiral knot 4i is practically 
suppressed and that the torus knot 5i becomes more fre- 
quent than the twist knot 52- This is an effect that was 
previously found for almost ideal configurations [1, Q and 
for polymers confined in spheres [HI ■ Here it is extended 
to the case of random self-attracting polymers in the col- 
lapsed phase. 

In previous sections we have shown that the aver- 
age writhe of a given knot is insensitive to changes 
in N. It thus makes sense to use the same windows 
[H^min, W'Vmax] for any comparison between results for 
different TV's. On the other hand, we have seen that the 
variance ctjv(t) 2 ~ N indicating that writhe distributions 
for contiguous (in complexity) knots will have an overlap 
that becomes more and more important as N increases. 
Hence, a windowing procedure should be less efficient in 
picking up a given subset of knots. This is indeed the 
case, for compact polygons with N = 1400 (filled his- 
tograms in Fig. [T3|) . For N — 1400 the values of the 
writhe considered do not show the writhe bias effect on 
the knot spectrum because the writhe distributions of the 
knots considered are not well separated. 

In the swollen phase, the conditional distributions 
P(t\ Wr) of different knots can be of comparable mag- 
nitude in a window of values [PFr m i n , Wr max ] only for 
very large N, because knots are rare for short chains. 
For long chains, P(t\ Wr) is essentially a broad Gaus- 
sian with average much smaller than its standard devia- 
tion, and of course with normalization ~ P(t). Hence, 
in a window [ Wr m { n , Wr max ] one just normally finds that 
P(r|[W r r m i n , Wr max ]) ~ P(t), that is, no sensible depen- 
dence on the writhe interval is present. 



IV. CONCLUSIONS 

We have investigated with Monte Carlo methods the 
interplay between writhe and knotting in ring polymers 
in both good and poor solvent conditions. To model good 
solvent conditions we used self-avoiding polygons on the 
simple cubic lattice and we incorporated an attractive in- 
teraction between neighboring pairs of vertices to model 
poor solvent conditions, leading to compact configura- 
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FIG. 13: Frequency of some knots, regardless of the writhe 
(upper panel), for compact configurations with 10 < Wr < 
20 (central panel), and 20 < Wr < 30 (lower panel), both 
for N = 600 (empty histogram) and for N = 1400 (dense 
histogram) . 



tions. 

For polygons with fixed knot type we computed the 
writhe distribution. The mean writhe is zero for achiral 
knots while for chiral knots it is non-zero but is insensitive 
to the length of the polygon. This has been observed 
previously for the good solvent case but we observe the 
same effect for poor solvents. Moreover, for each knot 
type, both mean values are consistent with that of ideal 
knots. 

In the swollen case, the width of the writhe distribu- 
tion increases as the square root of the length of the 
chain. Within errors, this is true also in the collapsed 
phase. On the other hand, the mean spread of the writhe 
computed without distinguishing knots scales differently 
(with exponent ~ 0.6) for compact polymers. The rea- 
son is that, in this regime, the population of knots grows 
quickly with the chain length, including complex knots 
with large mean writhe values. Therefore, the swollen 
phase is characterized by the geometrical spread of con- 
figurations, while the statistics of writhe in the compact 
regime is governed by the topological contribution of the 
mean writhe of knots. 

Finally, we have observed that picking up a specific 
knot by constraining the writhe in a given window is 
not possible if the chain length is much longer than the 
persistence length of the polymer. 

It is well known that the mobility of a circular DNA 
molecule in gel electrophoresis is sensitive to the writhe 
of the molecule. If the circular DNA is prepared with a 
particular knot type then this type of experiment could 
be used, in principle, to measure the mean writhe and 
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determine how the writhe is related to the knot type. 
By changing the ionic strength of the solution the typi- 
cal dimensions of the molecule could be changed, in an 
analogous way to changing solvent quality. 
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